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We present the first general search for new heavy particles, X, which decay via X —* WZ° — > 
ev + jj as a function of Mx and T(X) in pp collisions at ^Js =1.8 TeV. No evidence is found for 
production of X in 110 pb~ x of data collected by the Collider Detector at Fermilab. General cross 



section limits are set at the 95% C.L. as a function of mass and width of the new particle. The 
results are further interpreted as mass limits on the production of new heavy charged vector bosons 
which decay via W' — > WZ° in an extended gauge model as a function of the width, F(W'), and 
mixing factor between the W' and the Standard Model W bosons. 
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The Standard Model (SM) of particle physics is widely 
believed to be incomplete. Because alternative models 
are numerous and varied, it is advantageous to search 
for new physics using methods that are not specific to 
a single model, but which retain the most compelling 
aspects of currently favored scenarios [1] . 

A number of theories, including extended gauge mod- 
els, non-linear realizations of electroweak theory, a 
strongly interacting Higgs, and Technicolor, all pre- 
dict new high mass particles X which decay via X — ► 
WZ° [2-4]. A general X — > WZ° search can address 
all of these models, as well as new theories that may be 
proposed in the future. While typical searches in pp col- 
lisions, such as for a Technicolor pt — > WZ° [4], consider 
the narrow resonance case, T(X) <C Mj, there are good 
reasons to consider a general search which looks for X 
both as a function of mass and the width, even for large 
widths. For instance, a new heavy charged vector boson, 
W, has a width which can vary greatly as it depends on 
a mixing factor, £, between the W and the W [5]. 

In this Letter, we present the first general search for 
X -> WZ° production as a function of Mx and T(X) 
using pp collisions at \/s = 1.8 TeV using the Col- 
lider Detector at Fermilab (CDF). We use the final state 
WZ° — > evjj as it has the experimental advantages of a 
large branching ratio for Z° — > jj and a striking signa- 
ture of W -> Iv. The data, taken during the 1992-1995 
Tcvatron collider run, correspond to an integrated lumi- 
nosity of 110 pb _1 . Detailed descriptions of the detector 
can be found elsewhere [6] . The portions of the detector 
relevant to this search are: (i) a time projection cham- 
ber for vertex finding, (ii) a drift chamber immersed in 
a 1.4 T solenoidal magnetic field for tracking charged 
particles in the range \r\\ < 1.1 [7], and (iii) electromag- 
netic and hadronic calorimeters covering the pseudora- 
pidity range \rj\ < 4.2. An electron is identified as a 
narrow shower in the electromagnetic calorimeter that 
is matched in position with a track in the drift cham- 
ber. Jets are reconstructed as clusters of energy in the 
calorimet er using a fixed -cone algorithm with cone size 
AR = y/(A(j)) 2 + (At]) 2 = 0.4. The presence of neutri- 
nos is inferred from the momentum imbalance, ^t, in 
the transverse plane as measured in the calorimeters. 

Candidate events are selected online using a three-level 
trigger system [6] to identify W — > ev decays based on the 
requirement of an electron candidate with Et > 22 GeV, 
1 77 1 < 1.1 and a matching drift chamber track, and $t 
> 22 GeV. Several backup trigger paths, imposing for 
example electron Et > 25 GeV with no track require- 
ment and > 25 GeV, combine to make the trigger 
inefficiency for X — > WZ° production negligible in the 
final W — > ev sample. To search for resonant WZ° pro- 
duction, and to reduce standard model backgrounds, we 
raise the Et and $t thresholds and require two jets to 
be present. The final event selection requires an isolated 
electron [8] with E T > 30 GeV, $ T > 30 GeV, and two 
jets with Et > 50 GeV and 20 GeV respectively, each 
with \rj\ < 2.0. To reduce instrumental backgrounds, we 



restrict electrons to be in the fiducial region of the de- 
tector [9], and reject events in which significant hadron 
calorimeter energy is deposited out of time with the pp 
collision. A total of 512 events pass these requirements. 

The acceptance, Ax, for the process X — > WZ° — > 
evjj is defined as the number of events originating from 
X production and passing the final event selection, di- 
vided by the number of events in which X — > WZ°, 
W — ► ev; the Z° is allowed to have all decays. This 
definition allows non-quark decays of Z , such as Z° — ► 
t + t~ — ► jj, to contribute to the acceptance. To com- 
pute A x , we use the process W — > WZ° in the pythia 
Monte Carlo (MC) [10], followed by a parametric simu- 
lation of the CDF detector. We simulate W -> WZ° 
production for a variety of widths, T(W'). For narrow 
resonances, T(W') -C M\y>, the acceptance rises from 7% 
at M w , = 200 GcV/c 2 to 31% at M w , = 600 GeV/c 2 . 
For a given mass, the acceptance falls with increasing 
particle width. 

New production would show up as a resonance (peak) 
in both the dijet mass (M<jij e t = Mz») and the IF+dijet 
mass (Mv^+dijet = Mx)- Since a signal would appear as 
a clustered excess of events above the background spec- 
trum, we search by analyzing the shape of the data in 
the dijet vs. IF+dijet mass plane. Invariant masses 
are calculated using the measured energies and direc- 
tions of the electron, $t, and the two jets. To form 
the M^+dijet mass we fix the mass of the electron+i^T 
system to be equal to the W boson mass, which restricts 
the neutrino's unmeasured longitudinal momentum, p v z , 
to at most two possible values. When there are two so- 
lutions, we choose the one that yields a lower IF+dijct 
invariant mass. When there is no solution, we fix p v z 
such that the reconstructed W mass equals the trans- 
verse mass: M w = M T = 2p T p T [l - cos(</> e - <p v )} 1/2 . 
For T(X) <C Mx, MC studies show that on average these 
choices correctly reproduce the Z° and X masses with 
15% resolution and no significant bias. For a given T(X) 
the VF+dijet mass distribution is given by this mass res- 
olution and the intrinsic particle width. 

The primary background to this search is SM 
IF+jets production with W — > ev. To esti- 
mate this background, we use the vecbos MC [11] 
with Q 2 =<P T partons > 2 +M2 / , mrsdo' structure func- 
tions [12], herwig parton fragmentation [13], and the 
detector simulation. The W — > tv — > evvv background 
is similarly estimated but with TAUOLA MC [14] used to 
model the decay t — ► evv. We use vecbos to model the 
kinematics of the events, but use the data for an over- 
all normalization. Other backgrounds which produce the 
evjj final state include SM production of ti, W + W~ , tb, 
WZ°, Z°(-> e+e~) + jets, Z°(-* r+r") +jcts, and mul- 
tijet fakes. The ti, W+W~, tb, and WZ° backgrounds 
directly produce evjj events. Each is estimated using 
PYTHIA and the detector simulation, and is normalized 
to the measured or theoretical cross sections [15]. Wc 
estimate that there are 45 ± 14 ti events, 9 ± 3 W + W~ 



events, 3.0 ± 0.9 tb events and 1.6 ± 0.5 WZ° events in 
the data. The e+e~) + jets and Z°(-> t + t~) + 

jets events can fake the e$xjj signature if an electron 
> 2 jet event is lost, faking the 



from a Z' 



0l 



neutrino signature, or if in a Z° + 1 jet event an energy 
mismeasurement gives fake $t and an electron or tau is 
misidentified as a jet. We estimate the Z°+jets back- 
grounds using a combination of the PYTHIA and vecbos 
MC programs and the detector simulation, and normalize 
to the measured number of Z° + 1 jet data events in the 
Z° — > e + e~ channel. We estimate that there are 36 ± 5 
Z°(-> e+e~) + jets events and 1.6 ± 0.6 r+r") + 

jets events in the data. QCD multijet events can fake the 
evjj signature if a jet is misidentified as an electron and 
an energy mismeasurement in the calorimeter causes 
We estimate this background from the data in a manner 
similar to that used in Ref. [16], and predict that 27 ± 
3 QCD multijet events remain in the final sample. The 
contribution from all processes other than W+jets is 123 
± 16 events. 

We use a binned likelihood fit in the dijet vs. W+dijet 
mass plane (20 GeV/c 2 x 20 GeV/c 2 bins) to search 
for resonant WZ° production. All backgrounds except 
W+jets are normalized absolutely. The normalization of 
the W+jets background in the fit is fixed such that the 
sum of the signal and all backgrounds equals the number 
of events observed in the data. The relative magnitude 
of the signal and the PF+jets background is the only free 
parameter in the fit [17]. The VF+dijet mass spectrum 
for the data and background is shown in Fig. 1 for events 
with the dijet mass around Mz" and in the regions out- 
side a 25 GeV/c 2 mass window, with the expected dis- 
tributions plotted assuming no signal contribution. The 
results of the fit require no significant signal contribution 
and there is no evidence of resonant WZ° production for 
any mass or width for the acceptance model. 

To set general limits on the process pp — > X — > WZ°, 
we take X to be a W in an extended gauge model as it 
spans both the Mx and T(X) parameter space. Follow- 
ing the prescription of Ref. [5] (no additional fermions 
and the W and W vertex couplings, Wqq', Wiv and 
WWZ°, are identical), the production cross sections are 
uniquely determined as a function of mass, and the par- 
tial width of the W, T(W -> WZ°), is determined by a 
mixing factor, labeled £, which describes the amount of 
mixing between the W and the W . While this makes T a 
free parameter in the theory, we quote results in two spe- 
cific cases. The full mixing case, or reference model [5], is 
where the new particle X couples in the same way as the 
SMW{£= 1) and gives T{W -> WZ°) cx M^r, yielding 
a large branching fraction into WZ° , and widths compa- 
rable to the mass for Mw> w 425 GeV/c 2 . A second 
special case is £ — (4P^) 2 as in extended gauge models 
which restore left-right symmetry to the weak force and 
predict an effective W'WZ vertex term [2]. In this case, 
T(W) < M w > for all masses. 

We set limits on cr(pp — > X) ■ B, where B = B(X — » 
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FIG. 1. The VF+dijet mass spectra for the data and 
background with three different dijet mass requirements: 
Mdijot < 66 GeV/c 2 , 66 < M d ijet < 126 GeV/c 2 and 
Aidijct > 126 GeV/c 2 . A signal would appear as a resonance 
in the middle plot. The IF+jets background spectrum is nor- 
malized as described in the text so that there are the same 
number of events in the data as in the backgrounds. The up- 
per and lower plots show that the region outside the signal 
region is well modeled. 



WZ°) ■ B(W — > ev), using the fit technique described 
above and convoluting in systematic uncertainties, which 
depend on both mass and width, using the same methods 
as in Ref. [18]. The dominant source of uncertainty is the 
jet energy scale which would bias the measurement of the 
dijet and IT+dijet masses from the new particle X. The 
effect of such a bias is largest at lower mass, where in- 
creased background in the signal region can cause a large 
variation in the cross section limit. For example, the ef- 
fect is between 50% and 100% for the reference model and 
between 30% and 60% for £ = (^) 2 . Other notable 
sources of uncertainty are: uncertainty in the jet resolu- 
tion (between 15% and 30%), effect of the Q 2 scale on the 
IF+jets background shape (between 5% and 25%), choice 
of parton distribution functions (between 10% and 30%), 
uncertainty in W acceptance (between 5% and 30%), 
and MC modeling of initial and final state radiation (be- 
tween 5% and 15%). The total systematic uncertainty 
is found by adding the above sources in quadrature, and 
varies between 50-100% for the reference model and 40- 
75% for ^(^) 2 . 

The 95% C.L. upper limits on a ■ B for M x =200, 300, 
400 and 500 GeV/c 2 are shown in Fig. 2 as a function of 
the width. While these limits are not sensitive enough 
to set mass limits on pt — > WZ° production, they ex- 
clude a large region of W parameter space. Table I 
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FIG. 2. The 95% C.L. upper limits on a ■ B as a func- 
tion of the width. We use F(W // ) as it uniquely determines 
the W' — > WZ° branching ratio. Our results include the 
W' — > WZ° and W — > branching ratios. 



gives a summary of results for the r(W / ') <C M-^n ap- 
proximation using £ = (jjj^-) 2 - The results in Fig. 2 
can be interpreted as the first cross section limits as a 
function of W width, and Fig. 3 shows the first W 
exclusion region for £ vs. Mw where the theoreti- 
cal cross section exceeds the calculated 95% C.L. upper 
limit. Other direct searches for W 1 at the Tevatron in 
the W — > Iv and W' — > j'j channels have established a 
limit of M r > 786 GeV/c 2 [16,18,19], but only in the 
region of £ w which is complementary to region ex- 
cluded in Fig. 3, and are only valid for Y(W) <C My/>- 
A previous search for W' WZ° [20] sets cross section 
limits for M w , = 200,350 and 500 GeV/c 2 , but only for 
T{W')^M W ,. 

For a W in the reference model (£ = 1), we exclude 
the region 200 < M w , < 480 GeV/c 2 . For masses be- 
low 200 GeV/c 2 , the widths are small and the reference 
model is excluded at the 95% C.L. by the W' —* Iv re- 
sults [5,16]. Since the reference model is no longer valid 
for masses above 425 GeV/c 2 (T(W) becomes so large 
that perturbation theory is no longer valid [21]) the en- 
tire model is now excluded. These results are generally 
applicable to other new particles X with wide widths [3] . 

In conclusion, we have conducted a general search for 
new particles which decay via X — + WZ° in the eujj 
channel. We observe no evidence of resonant production 
and estimate production cross section limits as a function 
of mass and width. The results are further interpreted 
as mass limits on the production of new heavy charged 
vector bosons which decay via W — > WZ° in extended 
gauge models as a function of the width, r(W'), and 



TABLE I. Final results from the X -> WZ° search using 
110 pb _1 of data for T(X) < M x ■ Here we have modeled 
the new particle production with W' — > WZ° in an extended 



gauge model for the special case of £ = ( 



where £ 



is the mixing factor between the W and the SM W boson. 
Note that the 95% C.L. cross section upper limit from the fit 
is on X — > WZ° with W — * ev. In the denominator of the 
acceptance, Ax, we have allowed Z to have all decays. 



M x 


A x 


95% C.L. a ■ B Limit 


(GeV/c 2 ) 




(Pb) 


200 


0.07 


9.5 


300 


0.17 


4.5 


400 


0.24 


1.3 


500 


0.29 


0.7 


600 


0.31 


0.5 



mixing factor between the W and the W bosons. These 
are the first limits onl^ WZ° as a function of both 
mass and width, and are the only direct mass limits on 
W -» WZ° to date. 

We thank Jo Anne Hcwett, Tom Rizzo, and Jon Rosner 
for useful discussions and Jose Benlloch, Dan Hennessy, 
Marcus Hohlmann, Sacha Kopp, the Fermilab staff and 
the technical staffs of the participating institutions for 
their vital contributions. This work was supported by the 
U.S. Department of Energy and National Science Foun- 
dation; the Italian Istituto Nazionale di Fisica Nucleare; 
the Ministry of Education, Science, Sports and Culture 
of Japan; the Natural Sciences and Engineering Research 
Council of Canada; the National Science Council of the 
Republic of China; the Swiss National Science Founda- 
tion; the A. P. Sloan Foundation; the Bundcsministerium 
fuer Bildung und Forschung, Germany; the Korea Sci- 
ence and Engineering Foundation (KoSEF); the Korea 
Research Foundation; and the Comision Interministerial 
de Ciencia y Tecnologia, Spain. 



[1] Examples of such analyses include CDF Collaboration, 
F. Abe et al, Phys. Rev. Lett. 81, 1791 (1998); D0 Col- 
laboration, B. Abbott et al, Phys. Rev. Lett. 86, 3712 
(2001). 

[2] J. Pati and A. Salam, Phys. Rev. D 10, 275 (1974); 
R.N. Mohapatra and J. Pati, Phys. Rev. D 11, 566 
(1975); 11, 2558 (1975); G. Senjanovic and R.N. Moha- 
patra, Phys. Rev. D 12, 1502 (1975); P. Ramond, Ann. 
Rev. Nucl. Part. Sci. 33, 31 (1983). 

[3] P. Chiappetta and S. Narison, Phys. Lett. B 198, 421 

(1987) ; R. Casalbuoni et al, Nucl. Phys. B 310, 181 

(1988) . 

[4] K. Lane and E. Eichten, Phys. Lett. B 222, 274 (1989); 




200 250 300 350 400 450 500 550 600 
W 7 Mass (GeV/c 2 ) 



FIG. 3. The 95% C.L. excluded region in the £ vs. M w , 
plane, where £ is the mixing factor between the W and the 
SM W boson. While the branching fraction goes up with 
increasing values of £ and mass, the acceptance goes down as 
£ and the width increase. This causes the "nose" effect in the 
exclusion region. The largest mass exclusion occurs for £ = 
0.3, where we exclude M W i < 560 GeV/c 2 . 



D 50, 6734 (1994). 

[13] We have interfaced the subsequent parton shower evolu- 
tion and hadronization by the interface using the routines 
of G. Marchesini and B.R. Webber, Nucl. Phys. B 310, 
461 (1988); G. Marchesini et al, Comput. Phys. Com- 
mun. 67, 465 (1992). 

[14] S. Jadach, Z. Was, R. Decker, and J.H. Kuhn, tauola 
version 2.4, Comput. Phys. Commun. 76, 361 (1993). 

[15] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 73, 
2667 (1994); S.S.D. Willenbrock and D.A. Dicus, Phys. 
Rev. D 34, 155 (1986); J. Ohnemus, Phys. Rev. D 44, 
1403 (1991); 44, 3477 (1991). 

[16] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 74, 
2900 (1995). 

[17] We note that the assumption of fixing the normalization 
of the non-VF backgrounds has between a l%-5% effect 
on the limit. This variation has been incorporated into 
the overall systematic uncertainty. 

[18] CDF Collaboration, F. Abe et al, Phys. Rev. Lett. 74, 
3538 (1995). 

[19] CDF Collaboration, T. Affolder et al, submitted to 
Phys. Rev. Lett., hep-ex/0107008; D0 Collaboration, 
S. Abachi et al, Phys. Lett. B 358, 405 (1995); Phys. 
Rev. Lett. 76, 3271 (1996). 

[20] D0 Collaboration, V. M. Abazov et al, Phys. Rev. Lett. 
87, 231801 (2001). 

[21] J. Gunion, H. Haber, G. Kane, and S. Dawson, The Higgs 
Hunter's Guide, Frontiers in Physics series (1989). 



K. Lane, Phys. Rev. D 60, 075007 (1999); also see 
K. Lane, hep-ph/0007304 (2000). 
[5] G. Altarelli, B. Mele, and M. Ruiz-Altaba, Z. Phys. C 

45, 109 (1989). 

[6] CDF Collaboration, F. Abe et al, Nucl. Instrum. Meth- 
ods 271, 387 (1988). 

[7] We use cylindrical coordinates where positive z points 
along the proton beam and is zero at the center of 
the detector. The pseudorapidity, r\, is defined as r\ = 
— ln[tan(#/2)], where 9 is the polar angle with respect to 
the proton beam direction and <f> is the azimuthal angle. 
The transverse energy is defined as Et = Esm6, where 
E is measured in the calorimeter. 

[8] We require that the electron candidate pass identification 
and isolation requirements. The scalar sum of the pr of 
all tracks in the tracking chamber within a cone of Ai? = 
0.25 surrounding, but not including, the electron be less 
than 5 GeV/c, and that the Et in a cone of Ai? = 0.4 
around, but not including, the electron candidate be less 
than 10% of the electron Et- 

[9] The fiducial region is 0.05 < \rj\ < 1.05 and away from the 
edges of the calorimeter; for a more complete discussion 
of the fiducial region, see section 2.c. of F. Abe et al, 
Phys. Rev. D 52, 2624 (1995). 
[10] H. Bengtsson and T. Sjostrand, Comput. Phys. Commun. 

46, 43 (1987). 

[11] FA. Berends, W.T. Giele, H. Kuijf, and B. Tausk, Nucl. 

Phys. B 357, 32 (1991). 
[12] A.D. Martin, W.J. Stirling, and R.G. Roberts Phys. Rev. 



